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Hubbard chain, and a model in higher dimensions which exhibits the Mott 
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between different species. For the Bethe ansatz solvable models, the low tem- 
perature properties are well described by the excluson description if the degen- 
eracies due to string solutions with complex rapidities are taken into account 
correctly. For the Hubbard chain with strong but finite coupling, charge-spin 
separation is shown for thermodynamics at low temperatures. Moreover, we 
present an exactly solvable model in arbitrary dimensions which, in addition 
to giving a perspective view of spin-charge separation, constitutes an explicit 
example of mutual exclusion statistics in more than two dimensions. 
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I. INTRODUCTION 



Elementary particles or excitations are usually classified either as boson or as fermion. In 
recent years, however, it has been recognized that particles with "fractional statistics" inter- 
mediate between boson and fermion can exist in two-dimensional Jll or in one-dimensional 
systems. In two dimensions, a type of fractional statistics can be defined on the basis 
of the phase factor, exp{i6) with 6 allowed to be arbitrary, associated with an exchange of 
identical particles. One has ^ = for bosons and 6 = ir for fermions. A particle obeying 
such fractional statistics (with ^ 7^ or vr) is called as "anyon" |Q. It is believed that 
the quasiparticles and the quasiholes in the fractional quantum Hall liquids are anyons 0. 
Anyons can exist only in two spatial dimensions due to the braid group structure associated 
with them [Q. 

Another aspect of quantum statistics involves state counting, or the exclusive nature of 
the particles. Any number of bosons can be in a single-particle quantum state. Therefore 
there is no exclusion between bosons. On the other hand, the exclusion is perfect for 
fermions in a sense that a single particle state can accommodate at most one fermion. This 
aspect of quantum statistics can be generalized, as noticed by Haldane who proposed a 
definite generalization of Pauli principle such that one can consider particles with non-perfect 
exclusion. He pointed out that a spinon in one-dimensional long-range interacting quantum 
spin chain can exclude, on average, half of other spinon in occupying a single particle state. 
We shall call such a generalization as "exclusion statistics" and a particle obeying it as 
an "excluson". In contrast to usual bosons and fermions, the general concept of exclusons 
allows mutual statistics. Namely, there may exist statistical interactions or mutual exclusion 
between different species of particles. Haldane has recognized that quasiparticles in the 
fractional quantum Hall fluids are exclusons with mutual statistics between quasi-electrons 
and quasi-holes. 

Thus the concepts of fractional anyon statistics and exclusion statistics constitute gener- 
alizations of two different aspects (exchange phase and exclusion) of usual quantum statistics. 
An essential difference between the two concepts is that anyons can exist only in two spatial 
dimensions, while in principle exclusons may exist in any dimensions. 

Recently, one of us introduced the concept of generalized ideal gas of exclusons (see 
Sec. ^ for definition), and showed that its thermodynamic properties can be easily under- 
stood through a statistical distribution that interpolates between bosons and fermions. Later 
Bernard and Wu ^ have shown that the Bethe ansatz solvable models in one dimension 
can be described as an ideal (or non-interacting) excluson gas. This exemplifies that in 
certain circumstances particle-particle interactions can be totally absorbed by the statistical 
interactions. (A possible relation to a conformal field theory was discussed by Fukui and 
Kawakami |jl3|-) Thus the concept of exclusion statistics may become a powerful tool in 
understanding certain interacting many-body problems. For example, recently it has been 
shown P] that the essential features of low-temperature physics of Luttinger liquids in one 
dimension can be approximately described by a system of noninteracting exclusons. This 
may provide a new approach to interacting many-body systems. 



The examples considered by Bernard and Wu are the repulsive 5-function boson gas |10 



and the Calogero-Sutherland model |Tl|. Both of them contain only single species. The 
present paper is a follow-up to study statistical interactions or mutual statistics in exactly 
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solvable models with internal quantum numbers. We will discuss one-dimensional Bethe 
ansatz solvable models such as the isotropic (XXX) and anisotropic (XXZ) Heisenberg chain, 
and the Hubbard chain. In addition, we will consider an exactly solvable model in arbitrary 



dimensions proposed by two of us which has the Mott metal-insulator transition. A 
common feature in all these models is that there exists mutual exclusion either between 
particles of different rapidities or between different species. 

In this paper we are going to address the following two questions. First, many Bethe 
ansatz solvable models with internal degrees of freedom allow solutions of the Bethe ansatz 
equations with complex rapidities, while a naive analysis of Yang- Yang thermodynamic 
Bethe ansatz in terms of exclusons deals only with excitations which correspond to solutions 
with real rapidities. The question is whether these excitations are sufficient to give correct 



thermodynamic properties at low temperatures or not. Previously in ref. [0 the same topic 
of thermodynamics in terms of exclusons with multi-species is discussed without addressing 
this problem in the Bethe ansatz solvable models. For the Heisenberg spin chain, we ex- 
plicitly show that at the isotropic antiferromagnetic point (with the anisotropy parameter 
A = 1), the generalized ideal excluson gas of excitations with real rapidities does give cor- 
rect low-temperature behavior if a double degeneracy of the excitations induced by complex 
rapidities is taken into account. The second question deals with the excluson description for 
the physically interesting phenomenon of charge-spin separation. Both the Hubbard chain 
and the exactly solvable model in higher dimensions exhibit this phenomenon under certain 
conditions. We will show that when this happens, indeed the two models can be described 
by two species (spin and charge) of exclusonic excitations with nontrivial mutual statistics. 
Particularly the latter model provides an example of (mutual) exclusion statistics in more 
than two dimensions. We note that our treatment of the models differs from that presented 
in ref. M . 



II. EXCLUSON DESCRIPTION 

We consider a system with a total number = Nj' of particles or quasiparticles, where 
Nj' is the number of particles of species fi with a set of good quantum numbers, collectively 
denoted by j, specifying the states. Following [|^, we assume that the total number of states 
with {N^} is 

^=u \:,)iL.J .1 > (2-1) 



Nt\ Df{{m}) - 1 



where ({A^J}) is the number of available single particle states (counted as bosons), which 
by definition is given by 

^r({iv;})+E^r^; = G'r, (2.2) 

with statistical interactions gij'^ . Here Gf is the number of available single particle states 
when there is no particle in the system. Namely, Gf = ({0}). The derivative of (|2.2| ) is 

dD^iW"^}) 
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agreeing with the original definition for ^^statistical interactions" proposed by Haldane 0. 
When the pair {i, fi) differs from {j, u), we call g^^^ mutual statistics between particles labelled 
by {i, n) and those by (j, v). 

We further assume that the total energy of the system with {iVj^} particles is always 
simply given by 



(2.4) 



with constant e^. Equations ( |2.1| ), ( |2.2| ) and ( |2.4| ) define the generalized ideal gas of exclusons 
0. It is known that ( p.4|) is not satisfied for free anyons |T^. One of us has derived 
statistical distribution for generalized ideal gas and the thermodynamics following from it. 
The equilibrium statistical distribution for {A^f } is determined by 



where satisfy the equations 

(1 + <) n (y 



(2.5) 



exp 



T 



(2.6) 



where is the chemical potential for particles of species /i. The thermodynamic potential 
is given by 



VL = -TlogZ 
= -T^Gf log 



11,1 



C"^ -I- — V n'^^ /V 



(2.7) 
(2.8) 



where Z is the grand partition function. 



III. THERMODYNAMIC BETHE ANSATZ 

A large class of models which can be interpreted as generalized ideal gas is the Bethe ansatz 
solvable models in one dimension. It is well known that the generalized ideal gas actually 
represents a system of interacting particles. The Bethe ansatz method can be applied to 
systems with the special property that all the scattering amplitudes for many quasiparticles 
can be written in terms of two-body scattering amplitudes ||TB|. As a result, the eigenvalue 
problem of a Bethe ansatz solvable system is reduced to solving the so-called Bethe ansatz 
equations [|T^]. 

In general, the Bethe ansatz equations may be rewritten as 

Lp,{X^) = 2nl^ + Y: e,,{\l yj, (3.1) 

with 6fj,y{\1,\^ being the the two-body scattering phase shift between the quasiparticles 
with rapidities and X^. Here L is the size of the system, p^{\^) is a function of the 
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rapidity (or the quasi-momentum) , and {Ie}e,ii is a set of integers or half-odd integers. 
We assume that the integers {I^} satisfy > I^. We restrict ourselves to the solutions 
of Eq. ( |3.1| ) with rapidities all real. It is known that solutions with complex rapidities 
do exist for many of Bethe ansatz solvable systems. But we have not been able to deal with 
the problem of state counting for such solutions. Below we will address, in several concrete 
examples, the question of how important are the contributions of solutions with complex 
rapidities to low-temperature properties of the system. 

Subtracting the Bethe ansatz equations for adjacent rapidities A^ gives 



p.(A^+i) - pM)] = Mi^+i - If) + E [oAK+i, a;;) - e,,{\i Kn)] ■ (3.2) 



Following 1|18|,0 , we introduce a particle density 



L(A^ 



and a hole density 



^+1 



(3.3) 



(3.4) 



where M^^ > is the number of the holes in the branches of the Bethe ansatz equations 
(p.l[), i.e., the number is given by 

I^^, -1^ = 1 + M^. (3.5) 

In terms of these densities and we can rewrite the Bethe ansatz equations ( p.2|) as 

^p',{K)AXf = Lp,{Xf)AXf + Lp;(Ar) AAf + KiK, A,^) AAf] Lp.(Ap AA,^ (3.6) 



where we denote the derivative of a function / by /'. Here Lp^(Af)AAf and Lp|^(Af)AAf 
are, respectively, the number of the quasiparticles and the number of the quasihole in the 
interval AAf . The total number of possible choices of states with densities {p^} and {p^} is 
given by 



TBA 



n 



{L[p,(An+p;:(An]AAf}! 



[Lp,(Ar)AAr]! Lp;j(Ar)AAr 



(3.7) 



This number PVtba of states, obtained from the Bethe ansatz equations ( |3.6| ), is of the 
exactly the same form as Eq. ( p.lD for the number W of states in a generalized ideal gas of 
exclusons with the statistical interactions given by Eq. (|2.2|). In fact, we get the equivalence 
by setting 



L 



= 7^K(Ar)AAr 



(3.8) 



iVf = Lp,(Ar)AAf , D^{{m}) = Lp%X>t)AX>t 



(3.9) 
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and 

= 6^,6,, + ^^^;.(Ar, AJ)AAr. (3.10) 
In the thermodynamic hmit, the Bethe ansatz equations ( p.6|) become 

^P,W = PJ(A) + E / 9'''{K X')pA^')dX' (3.11) 
with the statistical interactions 

g'^'^iX, A') = 6,JiX - A') + ^^^;,(A, A'). (3.12) 



The key point is that the Bethe ansatz equations ( p.llD are equivalent to the equations 



of the statistical interactions in the excluson formalism. 

Here we stress once more that in the above we have considered only solutions with real 
rapidities. The equation ( |3.12| ) for statistical interactions refers only to such excitations. But 
its validity is obviously independent of whether these solutions of Bethe ansatz equations are 
complete or not. Though right now we do not know how to do state counting for solutions 
with complex rapidities, we feel it plausible that a bit more complicated excluson picture 
may still apply. 



IV. THERMODYNAMIC FORMALISM AT FINITE TEMPERATURES 

The thermodynamic formalism for Bethe ansatz solvable models at finite temperatures 
was first developed by Yang and Yang for the cases without internal degrees of freedom. 
This formalism has been justified, e.g. , for a one-dimensional Bose gas with repulsive point 
interaction |jl9|]. In this section, we briefly review a straightforward generalization of the 
Yang- Yang thermodynamic formalism in general setting with internal degrees of freedom, 
and show that it is actually the same as generalized ideal gas with the identification (|3.8|)- 

In the thermodynamic limit, the number of particles of species fi per volume is given by 



L 



dX p^(A). 



(4.1) 



From (|3.7] ) , we have the entropy S = log Wtba as 



S 



T = Y. [pm(A) + piiX)] log [p,(A) + p;(A)] - p,(A) logp,(A) - p;(A) logp;(A)} 



(4.2) 



We assume that the total energy per volume is expressed as 

E _ 



f = E/^Ap,(A)eJ(A) 



(4.3) 
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in terms of an energy density e^. As a rule, this assumption is always satisfied in Bethe 
ansatz solvable models. 

The thermodynamic potential Q at equilibrium with temperature T can be evaluated by 
minimizing 

n = E-J2a^N^-TS (4.4) 

with respect to the variation of the densities and p^. Here is the chemical potential. 
As a result, we have [|^,|[ 

^ = E / ^^(A) log [l + w-^\X)] , (4.5) 

where the functions are determined by the equations 

log[l + ^.(A)] - E / dX'g'^^iX', A) log [l + w-\X')] = fMz^. (4.6) 

The particle densities at equilibrium are determined by the Bethe ansatz equations 

^P',W = P,Ww,{X) + \')Pu{X')dX'. (4.7) 

The hole densities are given by p^ = p^w^. 

For the thermodynamics reviewed here, the following cautious remark is in order. We 



have followed the Bethe ansatz method to reduce the energy eigenvalue problem in 



the models to solving the so-called Bethe ansatz equations. However, it has not yet been 
generally proved that all the solutions to the Bethe ansatz equations provide the complete 



set of energy eigenstates |20]. In particular, it has been known in a number of cases that 
there are solutions with complex rapidities to the Bethe ansatz equations in addition to 
those with real rapidities. Though exact thermodynamics must deal with all eigenstates in 
a complete basis, low-temperature properties of a system, which are main focus of interests 
in many situations and in the present paper, might involve only a set of solutions that are 
not necessarily complete. It has been shown that the ground state at absolute zero always 
corresponds to a solution with all rapidities real. A sensible question is thus whether the 
solutions with real rapidities are enough for accounting for low-temperature properties of the 
system. A way of investigating this problem is to apply the thermodynamics reviewed above 
to some cases where the exact thermodynamics has been studied by ways that avoid the 
completeness assumption, and then compare the results to the exact ones. One well-known 
example of such cases is the Heisenberg spin chain, which we are going to study in next 
section. 



V. ISOTROPIC AND ANISOTROPIC HEISENBERG SPIN CHAINS 

In this section, we apply the formalism reviewed in the last section to the isotropic and 
anisotropic quantum spin-1/2 chains. On one hand, for the spin-1/2 isotropic Heisenberg 



7 



antiferromagnetic chain at low temperatures, we are able to explicitly demonstrate the 
magnon excitations (des Cloizeaux-Pearson-Faddeev-Takhtajan mode ||2l|,2^) in the exclu- 
son description. On the other hand, we will show that for the isotropic antiferromagnetic 
chain, the low temperature behavior obtained from a generalized ideal gas of the doubly de- 
generated or spin- 1/2 excitations with real rapidities does agree with the known exact results 
23| , |2^ which are obtained in other ways without invoking the completeness assumption of 



string solutions . Here we stress that for getting the correct results, it is necessary to take 
into account a double degeneracy of the excitations induced by complex rapidities. Strictly 
speaking, the solutions of the Bethe-Ansatz equations excitations with only real rapidities 
generally do not properly account for the low-temperature behavior, so other solutions with 
complex rapidities have to be included. This problem has not been addressed in the previous 
treatment [O] on mutli-component systems in one dimension. 



A. Spin-1/2 Isotropic (XXX) Heisenberg Chain 

The Hamiltonian of the spin-1/2 XXZ chain is given by 

n = i: [SJS^^, + S]Sl, + AiS]S]^, - 1/4); 



(5.1) 



with the periodic boundary condition S^+i = Si, where Sj is the spin-1/2 operator at site 
j, and A is the anisotropy parameter. 

First consider the isotropic (A = 1) Heisenberg chain. As is well known, for this model, 
Bethe |]T^ first reduced the eigenvalue problem of the Hamiltonian to solving the so-called 
Bethe ansatz equations 



M 



2L tan-^(2A„) = 2ttI^ + 2 ^an"^ (A„ - Xp) , (a = 1, 2, ■ ■ ■ , M), 

/3=1 

where M is the number of down spins. The energy eigenvalues are given by 



AI 



E 



r, i + 4Ar 

Comparing ( p. 21) with (|3.1| ), we have 

p{Xa) = 2tan"^(2Ac,) 



and 



^(A„,A/3) = 2tan-^(A„-A/3). 
Further, from ( ^.12 ), we get the statistical interactions between magnons as 



TTi + {x-x'y 



(5.2) 

(5.3) 

(5.4) 
(5.5) 

(5.6) 



Yang and Yang proved |^ that the ground state at T = indeed corresponds to a solution 
with all rapidities real, so that it is the same as the ground state of a generalized ideal gas 
with above statistical interactions. 
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B. Spin-1/2 Isotropic Heisenberg Chain at Low Temperatures 



By substituting the results in Section |V A| into ( [4. 51 ) and ( |4.6| ), we get the thermodynamic 
potential fl as 



where e = logw is determined by the equation 



l + e" 



-e{X)/T 



e(A) 



J" POO 



d\'-_ 



1 



1 + 4A2 TT 7-00 1 + (A - A 
Further, from ( [4.7|) , we have 

12 1 '•"^ 



1 + e 



-e{X')/T 



(5.7) 



(5.8) 



p(A) 1 + 



rfA'p(A' 



7rl + 4A2 nJ-oo '1 + (A-A')2 

First consider the zero temperature limit. Then ( |5.8|) and { f).9\ j become 

2 I 1 



(5.9) 



and 



6(A) 



P(A) 



rfA'- 



1 + 4A2 TT 7-00 1 + (A - A') 



1 2 1 

TT 1 + 4A2 vr 



t;A'p(A' 



1 + (A-A')^ 



(5.10) 



(5.11) 



Since these equations are, respectively, linear with respect to the unknown functions e and 
p, one can easily obtain the solutions as 



and 



e(A) 



P(A) 



TT 



2 cosh ttA 



1 



2 cosh vrA 



(5.12) 



(5.13) 



By substituting ( |5.13| ) into the right-hand side of (|5.9|), i.e., by iteration, we get the 
expression 



P(A) 



2cosh7rAl + e<^)/^ 



(5.14) 



at low temperatures. 

To clarify the physical meaning of ( |5.12| ) and ( [5. 141 ), we introduce a new variable k G [0, it] 

by 



sinh ttA = — cot k. 
In terms of k, we can rewrite (|5.12|) and ( ^.14|) as 



(5.15) 
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e(A) = Sk 



vr , 



sin k\ 



(5.16) 



and 



p{\)dX 



1 



dk 



27rl + e^fc/^' 



Equation ( p.l6| ) implies that hole-like excitations have the dressed energy 
which is nothing but the des Cloizeaux-Pearson-Faddeev-Takhtajan mode [^Ii|22 
known 



(5.17) 

7r| sin fc|/2, 
It is 



2T|] that the momentum of this mode, q G (— vr, vr], is related to the variable k by 




< g < vr; 
-TT < g < 0. 



(5.18) 



This double degeneracy is induced by complex rapidities [^, whose net effect is to give 
spin-1/2 to the excitations |2^. A certain explanation of this degeneracy will be given in 



the case of the spin-1/2 XY chain below. Equation ( ^.17|) implies that in terms of the 
dressed energy ( |5.16| ), these hole-like excitations obey Fermi statistics at low temperatures. 
In fact, by combining ( |5.17|) with ( |4.2|) and p|*, = p^w^ = p^e"/^, we have the expression of 
the entropy at low temperatures as 



S_ 
L 



r dq {f{e,) log f{E,) + [1 - f{E,)] log[l - f{e,)]} 

ZtX J-tt 



with the Fermi distribution function 



(5.19) 



(5.20) 



This result agrees with the free fermion descriptions of the XXX chain |25|. The entropy 
( ^.19| ) gives the specific heat C = 2T/3 per volume at low temperatures. The resulting 
specific heat is just the known exact specific heat [^,0]. Although we have dropped (string) 
solutions with complex rapidities to the Bethe ansatz equations, we have been able to get 
the correct low temperature behavior by taking into account the double degeneracy of the 
magnon excitations induced by complex rapidities. Strictly speaking, the contributions to 
the specific heat from the string solutions cannot be ignored at low temperatures. But the 

of the 



total contribution |2^ appear to be just doubled, due to the spin-1/2 character 
magnon excitations. This implies that the low energy excitations as quasiparticle can be 
indeed described by the generalizad ideal gas of exclusons with real rapidities only. 



C. Spin-1/2 Anisotropic (XXZ and XY) Chains 

Next consider the spin-1/2 chains ( p.lj ) with anisotropy — 1 < A < 1. Then the Bethe ansatz 
equations are given by |^ 

M 
m=l 
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for the states with M down spins. Here 



and 



-2 tan"^ 



(cot T]) tanh [ ^ — ^ 



(5.22) 



(5.23) 



with A = cos 7], (0 < ?7 < tt) . The energy eigenvalues are given by 



M 



^ = -5][A-cosp(A,)]. 



(5.24) 



Comparing ( [5.21| ) with (p.l|), and using ( |3.12| ), we get the statistical interactions between 
magnons as 



(7(A,A') = 5(A-A') + 



1 



sin2r7 



2tx cosh(A — A') — cos 2rj 



(5.25) 



Again, the ground state at T = is a solution with all rapidities real ||2^, so that it coincides 
with the ground state of a generalized ideal gas with above statistical interactions. 
In the special case of the XY model, i.e., A = (77 = 7r/2), we have 



(5.26) 



This implies that quasiparticles (magnons) obey Fermi statistics. Of course, this is consistent 
with the well-known fact that the spin-1/2 XY chain can be transformed into a system of 
free spinless fermions in one dimension by using the Jordan- Wigner transformation p8| . 

Note that — (vr — rf) < p{\) < {n — rj) < > —00 < A < +00 from ( ^.22| ). In particular, 

we have -n/2 < p(A) <it/2 for the XY chain. Combining this with (|5.24|) , (^) and (^3D, 
we have the thermodynamic potential 



— = / (iplog {1 + exp[cosp/T]} 

L 2tT J~tt/2 



(5.27) 



for the XY chain. However this does not coincide with the known exact thermodynamic 
potential which is given by replacing the range of integral (|5.27|) to [— tt, vr]. One possible 
reason for this discrepancy is that we have not taken into account solutions with complex 
rapidities A^ to the Bethe ansatz equations (|5.21|) [1^. But, as is well known, there is no 
complex rapidity in the XY chain Where does this inconsistency come from ? The 

answer is the following. In the procedure to take the XY limit A — >■ 0, we missed a double 
degeneracy of the magnon excitations. Actually a more careful treatment taking into account 
complex rapidities shows this double degeneracy in the the XY limit P9 . A similar situation 



occurs also in the XXX chain at low temperatures [Q. Of course, a much simpler way to 
get the exact result of the XY chain is to use the variable p instead of the rapidity A. Then 
one can easily obtain the exact result. 
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VI. THE HUBBARD CHAIN 



In this section we turn to the Hubbard chain. Our intention is to present an exclusonic 
description for spin-charge separation, which is known to occur in this model with infinite 
couphng at zero temperature |^0[- We will demonstrate spin-charge separation under the 



following two broader conditions: (i) The on-site Coulomb energy U is large but finite; 
(ii) The temperature T is finite but sufficiently low compared to the order of the effective 
exchange 1/U. 

A. General Considerations 

The Hamiltonian of the Hubbard chain is given by 

L L 

^ = E E (S>Ci+i.- + Ci+i,<xCi,a) + UY^ c]^^9,tc]_|Cj,| (6.1) 

with on-site Coulomb energy U and periodic boundary conditions CL+i,a = ci^a (o^ =T)i), 
where ci „, respectively, the creation and annihilation operators for an electron with 

spin a at the site j. 



For the Hubbard chain, Lieb and Wu obtained the Bethe ansatz equations 

Lk'r = 27rI^ + Y.9(-smk':~yA, (z = 1, 2, ■ ■ ■ , iV) (6.2) 

/3=1 ' 

and 



-E^(-^^-^sinfc^^j =27r/^-^^^^^[^-^), (a = l,2,---,M) (6.3) 



with quasi-momenta {/cf}^^ and rapidities {\^a\^=\ (iV — M) electrons with up spin and 
M electrons with down spin. Here 

e[x) = -2iwi-^{2x). (6.4) 

The energy eigenvalues are given by 

N 

E = -2j2cosk''j. (6.5) 
i=i 

Comparing ( |6.2|) and (|6.3| ) with (|3.1|) , we get 



Pcik'r) = kt, (6.6) 

PsiK) = 0, (6.7) 
e,,{km) = Q, (6.8) 
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(6.9) 



^ss(A^, A^) 



and 



2 2 



MA:,fcj) = ^(A:-^sinA:j). 
Further, from ( |3.12D , we have 



(6.10) 



(6.11) 



(6.12) 



cos k 



TcU l + 4(2smA;/f/- A) 



2 ' 



(6.13) 



and 



r(A,A') = 5(A-A') + 



TT 



1 + (A-A02 



r (A, fc) 



1 



vr 1 + 4(A -2smk/Uy 



(6.14) 



(6.15) 



Since in the ground state all kj and A* are real the ground state is the same as 
that for a generalized ideal gas with two species with above statistical interactions. In ref. 
1 31 1 it has been shown that solutions with complex quasi- momenta and rapidities in the 
Hubbard model correspond to gapful excitations. Therefore at low temperatures, we can 
ignore thermal activations of such excitations and concentrate on solutions with only real 
quasi-momenta and rapidities. 

We note that our results for the Hubbard chain disagree with those in a previous paper 
on the same subject. The disagreement comes out from the incorrect procedure taken 
in Ref. |jl^ . In the second equation of Eq. (4) in Ref. 



the first and fourth terms cancel 
each other exactly. Namely, the fourth term in R.H.S., — J2^Li A/j^a^^.a^ = — Aq, cancels the 
first term. Therefore if one proceeds correctly, the first and fourth terms on the R.H.S. of 
the second equation of (5) must also cancel. For the same reason, the first and fourth terms 
on the R.H.S. of (9) should cancel as well. But surprisingly it is not the case in Ref. 
Two originally cancelling terms become eventually non-cancelling! 



B. Strong Coupling Expansion for the Hubbard Chain 

It is well known that spin and charge degrees of freedoms are separated in the strong coupling 
Hubbard chain at zero temperature PO] . To explore the possibility of a similar spin-charge 
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separation at finite temperatures, we first perform the 1/U expansion for tlie tliermodynamic 
potential of the Hubbard chain. 

To begin with, we will write down the explicit form of the thermodynamic potential Vl 
of the Hubbard chain. Substituting the results of the statistical interactions in Section |Vl7\| 
into the thermodynamic formulas ( [4.5|) and (|4.6| ) in Section fV|, we have the thermodynamic 
potential 

VL T 



dk log 



1 + e 



L 2tx 

for the Hubbard chain, where = \ogw^ are determined by the equations 

1 



2T roo 

tc{k) + — 

TT 



d\'- 



and 



+ 



. 1 + 4[A'- (2sinA;)/t/]2 

J" roo 



l0£ 



-2 cos k — ttr 



(6.16) 



(6.17) 



1 



TT J-oo 

4T 



dk 



1 + (A - A') 
cos k 



ttU J-n 1 +4[(2sinA;)/f/- A]2 



1 + e-^»(^')/^ 
log 



l + e" 



,(fc)/T 



(6.18) 



Now we proceed to the 1/U expansion. When we neglect higher orders than 1/U, the 
equations (|6.17|) and ( |6.18| ) become 



ec{k) = —2 cos k — Qc — 



2T 



TT 



1 



oo 1 + 4A2 



e.(A) 



2J 



+ 



J" roo 



1 + 4A2 vr J-oo 
with the effective exchange integral 

2T 



d\' 



1 + (A - A')' 



J = / dk cos k log 

7r(7 



l + e 



log 
log 

c(k)/T 



1 + e 



1 + e 



(6.19) 
(6.20) 

(6.21) 



of the spin-1/2 isotropic Heisenberg 



The equation ( |6.20[ ) for corresponds to Eq. 
chain, although es and ec couple to each other. 

Similarly we have the equations for the charge density pc and the spin density ps as 



Pc{.k) 



and 



1 + e 



ecik)/T 



^ + ^7 / d\ps{\) 



2n nU 
1 N 4 



4 



1 + 4A2 



1 



27r L 1 + 4A2 TT J-oo 



d\'ps{\') 



cos k 



1 



1 + (A - A')^ 



(6.22) 



(6.23) 



from (|4.7|) , where N is the number of electrons. Here we have neglected higher orders than 
1/U . Comparing ( |6.23| ) with ( p.9| ) of the spin-1/2 isotropic Heisenberg chain, we obtain that 
the number p^ of down spins per volume in the Hubbard chain is proportional to the number 
of the electrons per volume. 

Clearly, from ( p.l9| )- ([6.23| ), we conclude that the degrees of freedoms of spin and charge 
parts still couple to each other at finite temperatures. 
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C. Spin-Charge Separation in the Hubbard Chain at Low Temperatures 



In this section, we will demonstrate that the spin-charge separation occurs in the strong 
coupling Hubbard chain at low temperatures. 

Consider the case with low temperatures T ^ 1/U. That is to say, temperatures are 
sufficiently low compared to the effective exchange J ( |6.21| ) which behaves as 



J 



7tU Jd 



dk cos k[—ec{k)] 



(6.24) 



at low temperatures. Here = {k\ec{k) < 0}. Then we have = 0(1 /U) from (|6.20 
Note that ec can be rewritten as 



6e(A;) = ef(A;) + e 



.(1) 

-c 



with 



and 



-2 cos k — ar 



dX 



TT J-oo 



1 + 4A2 



log 



1 + e 



(6.25) 



(6.26) 



(6.27) 



Here e^^) = 0(1 /U) due to the above e, = 0{1/U). 

From this observation, we devide the integral of (|6.16|) into two parts as 



where 



T 















/ dk log 




Jd± 





(6.28) 



(6.29) 



and the ranges of integration are: = |A;|e^°^(A;) + e^^-*} > o| and 

D_ = {A;|e(°)(fc) + e« < o}. Since >> e« = 0{1/U) and T « 1/U, I± dOOl ) are 
given up to the order 1/U as 



T 



/ dk log 




Jd+ 





(6.30) 



and 



dklog 



1 + J 
1 + e' 



L ' 2tt 



l + e 



ef\k)/T 



(6.31) 
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where we have used the relation dk/ 271 = N/ L which is derived from (|6.22| ) at low 
temperatures and for large U. Thus we have 



-^--/ .Hog 



N2T dX 



TT 



oo 1 + 4A2 



log 



1 + e 



(6.32) 



for the thermodynamic potential Q ( p.l6| ) up to the order 1/U, where we have used ( |6.27D , 
and the thermodynamic potential Qc of the charge part is given by 



L 



T 
2^ 



dklog 



(6.33) 



Thus the thermodynamic potential Q of the Hubbard chain can be written as the sum 
of the charge and the spin parts 



where 



J" POO 

"2^ 



log 



1 + e 



(6.34) 



(6.35) 



which is identical to the thermodynamic potential ( |5.7|) of the spin-1/2 isotropic Heisenberg 
chain. Eq. ( |6.34D is an indication of charge-spin separation. 



VII. EXACTLY SOLVABLE MODEL IN HIGHER DIMENSIONS 

So far we have discussed the models in one dimension to which the Bethe ansatz approach is 
applicable. In two dimensions, it is known that the quasiparticles of the fractional quantum 
Hall liquid are anyons. They can also be considered to be exclusons. In this section we 
present an example of mutual exclusion between different species in an exactly solvable model 
in higher dimensions, that exhibits charge-spin separation under certain circumstances. This 
clearly shows that exclusion statistics is conceptually different from anyon statistics whose 
existence requires two (spatial) dimensions. 



Recently two of us [0 and Baskaran ||33| have proposed a model of interacting electrons 



that can be solved exactly in any dimensions. The Hamiltonian is 

U_ 

(i.j) ^ i,j,i,m 



H = - J2 4aCi,a + h.C. + J^J2 ^i+j,i+mcl^Cj,^cl^Cm,l - + (^/^o/i)4<,Q,<„ (7.1) 



where {i,j) represents nearest neighbors in d dimensions, and L"^ is the total number of 
lattice sites, fi the chemical potential, fiQ the magnetic moment and h the external magnetic 
field. 

This model is unrealistic in the sense that the interaction term with coefficient U is 
of infinite-ranged in real space and of strength independent of distance. (Note that the 
noninteracting electon models are even more unrealistic since they neglect the long-range 
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Coulomb.) It, however, has an attractive feature of being exactly solvable. In fact, it can 
be easily diagonalized for each k in momentum space. All the properties including the 
thermodynamic quantities were obtained in [0. Furthermore it is remarkable that this 
model exhibits a number of important features of the correlated electron problems in spite 
of its simplicity and unrealistic nature. It exhibits, for example, the Mott metal-insulator 



transition that was stressed by two of us |]T2[ and Continentino and Coutinho-Filho [p^j . 

The zero temperature phase diagram of this model in any dimensions is shown in Fig. ^ It 
has both the fixed-density and density-driven Mott transitions. These transition in general 
may be in different universal classes [^,^. The critical exponents of the two types of 
transitions are, however, the same in the present model and they seems to be in the same 
universality class |3^. In the zero temperature phase diagram Fig.l, the region OBC is 
a Mott insulator phase with half filled band. The rest is metallic phases. In the region 
OABC, a double occupancy is prohibited and, as we will show, the system can be described 
by an excluson picture. There is a Fermi surface of the excluson gas in the region OAB , 
and on the phase transition line OB is a quantum phase transition of the excluson gas. On 
the other hand, the excluson description breaks down on the phase transition line BC. The 
region outside OABC is a metallic phase which is described by the two species of fermions 
(spin up and spin down electrons) as it should be. 

Let us concentrate on the region OAB. Assume that U is large and T is low, so that U 
is much larger than both T and the band width [^]. Under these conditions, there is no 
activation of doubly occupied states, therefore there are only three states (0,1 and 2) for 
each momentum k. In the state there is no electron, the state 1 an electron with spin up, 
and the state 2 an electron with spin down. Let us denote the number of charges as Nc and 
the number of magnons (number of spin-down) as iV^. We regard and Ng as independent 
variables (spin-charge separation). By definition, the state has Nc = and A^^ = 0, the 
state 1 A^c = 1 and Ng = and the state 2 A'c = 1 and A^^ = 1 (See table |lD. Then from 



we easily derive 



I . Gg = 0, 



and 



'9cc 


9cs' 




" 1 0" 


.9sc 


9ss. 




.-1 1. 



(7.2) 



(7.3) 



It is easy to verify that the condition for the ideal excluson gas ( ^.4[ ) is satisfied, so that 
the system with doubly occupied states suppressed can be described as a generalized ideal 
gas with two (charge and magnon) species with the statistical interaction given by ( [7.31) . 
Note the nontrivial value —1 for mutual statistics Qsc'-, i-e. the presence of a charge can creat 
a magnon state, though there is no bare available single magnon state [Gg = 0) when there 
is no charge. It is straightforward to check that the thermodynamics of the generalized ideal 
excluson gas obtained from Eqs. (|2.6| ) and (|2.8| ) is identical to the result of Ref. |jl2|] in the 
low temperature limit. 

Indeed, in the present case, the species index fi = c,s, and the state index j is the 
momentum k in dimensional space. Equation (p^.8|) now takes the form 



d'^k 



log 



(k) 

1 - Y.u9fMunuik) 



(7.4) 
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where n^(/c) is the occupation number distribution function of the charge (/i = c) or spin 



(7.5) 



(/i = s) excitations in /c-space. From Eqs. (|2.5| ) and (|2.6|) , we have 
nc{k)[l + Wc{k)] = 1, ns{k)[l + Ws{k)] = ndk), 
where the statistics matrix ( [7.31) has been used , and Wc{k), Ws{k) satisfy 



Wc{k) 



l + Wsjk) 
Ws{k) 
Ws{k) 



,(e,(fc)-Mc)/T 



,(e4fc)-M,)/T 



(7.6) 
(7.7) 



where ec{k) = -2J2i=iCos{ka), e.^ 



measured relative to the energy of spin-up electrons) , and 
Thus 



2fiQh (the energy of spin excitation, which is actually 

(7.8) 



Wcik) 



,(€c(fc)-Mc)/T 



1 + (,-(es{k)-t^s)/T- 

Substituting Eqs. O, ( [TTD , and (^ into ([7^ ), we obtain 

d'^k 



(7.9) 



n 



-T 
-T 



(2.)^'^°S 
d'^k , 
(2^^°S 



1 + (^gMoVr ^ g-MoVr-)g(Mc-ec)/T 



g-(£i(A:)-/.)/T _^ g-(e2(fc)-M)/T 



(7.10) 



where ei{k) = tcik) — ^oh, and e2{k) = edk) + n^h are the energy of spin-up and spin- 
down electrons respectively. Equation (|7.10|) is nothing but the result of Ref. [1^ in the low 
temperature limit. (Remember that here we consider the case with large U and low T, so 
that doubly occupied states are suppressed.) 

Here we emphasize that the concept of spin-charge separation is crucial. The effects that 
spin and charge excitations are not actually independent of each other have been taken care 
of by the statistical interaction or mutual statistics between them in the present formulation. 
It seems to us that similar situations may happen in other strongly correlated systems that 
exhibit charge-spin separation in higher dimensions. We finally note that our treatment for 
the exactly solvable models in higher dimensions disagrees with that in a previous paper 
|M on the same subject. 



VIII. ACKNOWLEDGEMENTS 

This work was supported by Grant-in- Aid from the Ministry of Education, Science and 
Culture of Japan. The work of Y.S. Wu was also supported in part by U.S. NSF grant PHY- 
9309458. He is also grateful for warm hospitality of the Institute for Solid State Physics, 
University of Tokyo during his visits. 



18 



REFERENCES 



T electronic mail address: hatsugai@coral.t.u-tokyo. ac.jp 

electronic mail address: koma@riron.gakushuin.ac.jp 
[1] See, e.g., Fractional Statistics and Anyon Superconductivity, edited by F. Wilczek 

(World Scientific, Singapore, 1989) and references therein. 
[2] C. N. Yang and C. P. Yang, J. Math. Phys. 10, 1115 (1969). 
[3] F. D. M. Haldane, Phys. Rev. Lett. 66, 937 (1991). 
[4] F. Wilczeck, Phys. Rev. Lett. 48, 1144 (1982). 
[5] B. I. Halperin, Phys. Rev. Lett. 52, 1583 (1984). 
[6] Y. S. Wu, Phys. Rev. Lett. 52, 2103 (1984). 
[7] Y. S. Wu, Phys. Rev. Lett. 73, 922 (1994). 

[8] D. Bernard and Y. S. Wu, in Proc. of Sixth Nankai workshop (Tianjin, China; June 

1994), ed. by Mo-lin Ge and Y.S. Wu (World Scientific, 1995). 
[9] Y.S. Wu and Yuc Yu, Phys. Rev. Lett. 75 (1995) 890. 
[10] D. C. Mattis and E. H. Lieb, J. Math. Phys. 6, 304 (1965). 

[11] F. Calogero, J. Math. Phys. 10, 2191, 2197 (1969); B. Sutherland, J. Math. Phys. 12, 
251 (1971). 

[12] Y. Hatsugai and M. Kohmoto, J. Phys. Soc. Jpn. 61, 2056 (1992). 

[13] T. Fukui and N. Kawakami, Phys. Rev. B 51,5239 (1995). 
[14] K. Byczuk and J. Spalek, Phys. Rev. B 50, 11403 (1995). 
[15] Y. S. Wu, Phys. Rev. Lett. 53, 111 (1984). 

[16] See, e.g., B. Sutherland, in Exactly Solvable Problems in Condensed Matter and Rela- 
tivistic Field Theory, edited by B. S. Shastry, S. S. Jha and V. Singh, Lecture note in 
Physics. 242 (Springer, Berhn, 1985) p. 1. 

[17] H. A. Bethe, Z. Physik 71, 205 (1931). 

[18] E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 

[19] T. C. Dorlas, J. T. Lewis and J. V. Pule, Commun. Math. Phys. 124, 365 (1989); 

T. C. Dorlas, Commun. Math. Phys. 154, 347 (1993). 
[20] T. Koma and H. Ezawa, Prog. Theor. Phys. 78, 1009 (1987). 
[21] J. des Cloizeaux and J. J. Pearson, Phys. Rev. 128, 2131 (1962). 

[22] L.D. Faddeev and L. Takhtajan, Phys. Lett. 85A, 375 (1981); J. Sov. Math. 24, 241 
(1984). 

[23] M. Takahashi, Prog. Theor. Phys. 50, 1519 (1973); K. Nomura and M. Yamada, Phys. 
Rcv.B 43, 8217 (1991); J. Suzuki, Y. Akutsu and M. Wadati, J. Phys. Soc. Jpn. 50, 
2667 (1990); C. Destri, and H. J. dc Vega, Phys. Rev. Lett. 69, 2313 (1992). 

[24] H. M. Babujian, Nucl. Phys. B215[FS7], 317 (1983). 

[25] T. Yamada, Prog. Theor. Phys. 41, 880 (1969); L Affleck, Phys. Rev. Lett. 56, 746 
(1986). 

[26] T. Koma, Prog. Theor. Phys. 78, 1213 (1987); 81, 783 (1989). 
[27] C. N. Yang and C. P. Yang, Phys. Rev. 150, 321 (1966). 

[28] E. Lieb, T. Schuhz and D. Mattis, Ann. Phys. 16, 407 (1961); S. Katsura, Phys. Rev. 

127, 1508 (1962); 129, 2835 (1963). 
[29] M. Takahashi and M. Suzuki, Prog. Theor. Phys. 48, 2187 (1972). 
[30] M. Ogata and H. Shiba , Phys. Rev. B 41, 2326 (1990). 
[31] F. Woynarovich, J. Phys. C: Sohd State Phys. 15, (1982) 85, 97. 



19 



[32] E. H. Licb and F. Y. Wu, Phys. Rev. Lett. 20, 1445 (1968). 
[33] G. Baskaran, Int. J. Mod. Phys. Lett. B 5 643, (1991). 

[34] M. A. Continentino and M. D. Coutinho-Filho, Solid State Commun. 90, 619 (1994). 
[35] M. A. Continentino, Physics Report, 239, 179 (1994). 

[36] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. Fisher, Phys. Rev. B 40, 546 

(1989). 

[37] J. Spalek, and W. Wojcik, Phys. Rev. B 37, 1532 (1988). 



20 



FIGURES 

FIG. 1. Phase diagram of the exactly solvable model in higher dimensions at T = 0. 
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TABLES 



TABLE I. Electronic states and spin charge labels 



Electronic State 


Label 


Nc'. Charge 


Ng'. Magnon 














T 


1 


1 







2 


1 


]. 
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